Quasi phase matching for high order harmonic generation induced by the 

carrier-envelope phase. 
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We report a novel quasi-phase matching technique for high-order harmonic generation in low- 
density gases. Numerical simulations show that in few-optical cycle pulsed Bessel beams it is possible 
to control the pulse envelope and phase velocities which in turn allows to control the carrier-envelope 
phase during propagation. The resulting oscillations in the peak intensity allow to phase-match the 
high-order harmonic generation process with a nearly two decade enhancement in the XUV power 
spectrum. 
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Recent developments in ultrashort laser pulse technol- 
ogy have proved the viability and effectiveness of produc- 
ing Extreme Ultra- Violet (XUV) radiation by High Or- 
der Harmonic Generation (HHG) in low pressure gases. 
This process may be qualitatively understood through 
the three-step model P, by which the atomic poten- 
tial is strongly modified by the impinging laser radiation 
such that the electron may tunnel out of the potential 
well. Under the external influence of the strong driving 
laser field the electron wave-packet will be first driven 
away and then pulled back to the atom. The recombina- 
tion process will liberate the excess energy in the form of 
high energy XUV photons. Both the phase of the driving 
laser field and phase accumulated during propagation will 
finally determine how the XUV field coherently builds up 
so that in general, in order to achieve efficient HHG along 
a defined propagation direction phase matching require- 
ments must be fulfilled [1, [3j [B| ■ The propagation dis- 
tance over which this occurs is defined as the coherence 
length £c = 7r/|AA'|, where the total phase mismatch 
AK for HHG from the fundamental laser frequency to to 
the q'^ harmonic may be written as 



AK = A'disr 
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The contribution due to gas dispersion is /iTdisp = 
^fciasor — kq, whcrc fciascr and kq are the laser and Har- 
monic pulse wave-vectors, respectively. This term, which 
typically has positive sign, for low gas pressures (e.g. < 
50 mbar) can be neglected [1]. The plasma dispersion 
term may be written as i^^piasma = ^p(l ~ '7^)/29cw, 
where ujp is the plasma frequency, and gives a neg- 
ative contribution while the atomic dipole phase 
term -f^atomic, related to the quantum paths of the elec- 
trons involved in the HHG process, in general does not 
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have a fixed sign during propagation since it is pro- 
portional to the intensity gradient along the propaga- 
tion direction, dl/dz [1, Finally, the geometric term 
-K^geomotric dcpcuds On the shape of the driving laser pulse. 
For example the Gouy phase shift in a tightly focused 
Gaussian pump pulse will give a positive contribution 
^'geometric = 2((7— l)/b whcrc b is the confocal parameter 
whilst a pump pulse confined in the lowest order hollow- 
core waveguide mode will give a negative contribution 
^'^gcomotric = — 2c/2ti;a^, where un is the first zero of 
the Jq Bcssel function and a is the fiber radius 

Truncated Bessel beams, obtained by refocusing with 
a lens the output of a hollow core waveguide, have been 
proposed and the specific phase profile, related to the 
geometry of the driving pulse, has shown to lead to 
improved HHG brightness More recently Pulsed 

Bessel Beams (PBB) without successive refocusing, have 
also been proposed to optimize phase-matching in the 
XUV region [a, 113. PBBs have the pulse wave vec- 
tors aligned along a conical surface such that the result- 
ing interference pattern is characterized by a central in- 
tense peak with surrounding low-intensity rings llj . By 
changing the angle 7 of the wave-vector cone it is pos- 
sible to continuously tune the negative geometric term, 
-f^goomotric = " 'z('^/c)(l — C0S7) [9| and thus achieve per- 
fect phase-matching {AK ~ 0) for a wide range of oper- 
ating parameters. Due to the negative sign of this last 
contribution, phase-matching is possible only in the pres- 
ence of a positive material dispersion (e.g. at high gas 
pressures) that is larger than the negative plasma contri- 
bution. 

In the absence of the possibility to achieve perfect 
phase-matching the next best thing is so-called Quasi- 
Phase-Matching (QPM) by which the HHG process is 
periodically strongly reduced or eliminated at propaga- 
tion distances corresponding to odd multiples of the co- 
herence length ic- This fruitful technique has been im- 
plemented by a number of means some of which involve 
a static modulation of the gas medium while others are 
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FIG. 1: (in color online) (a) Scheme of pulsed Bessel beam 
geometry. 7 and 5 are the Bessel con angle and the pulse 
front tilt angle, respectively. 'z=0' indicates the position of 
the input condition for the numerical simulations, (b) shows 
the calculated peak intensity evolution on axis for the PBB 
for two different cone angles, 7=1, 1.5 deg. 



based on a more flexible and tunable counter-propagating 
beam configuration 

In this work we propose a novel technique to achieve 
QPM with few-cycle pulses by controlling the Carricr- 
Envelope-Phase (CEP) along propagation. This tech- 
nique is very general and may in principle be applied 
to any nonlinear intensity-dependent process although 
here we consider the specific case of HHG. The under- 
lying idea is very simple: when the pulse contains only a 
few optical cycles the instantaneous peak intensity /inst 
will become very sensitive to the CEP value, i.e. to the 
relative phase between the carrier- wave and the pulse en- 
velope. For CEP=0 (7r/2), /jnst will be maximum (min- 
imum). Since HHG is extremely sensitive to variations 
in the pulse intensity and phase, a periodic modulation 
of the instantaneous intensity and phase may lead to a 
quasi-phase-matched conversion process. 

This technique requires control over both the phase 
velocity, which determines how the carrier-wave propa- 
gates, and the group velocity, which determines how the 
pulse envelope propagates. Our proposal is to achieve 
this by use of PBBs. We first note that the PBB is just 
a particular case of a larger family of "conical" pulses for 
which the wave-vectors propagate along a cone with an- 
gle 7 and a pulse front that in general, may be tilted by 
an angle S with respect to the wave- vector. The carrier- 
wave velocity, w^, and the envelope peak velocity, Vp, 
along the propagation direction are related to the cone 
and tilt angles by 



W0 = Wi/(iCiC0S7) 

Vp = cos 6 / K' cos{^ + S) 



(2) 
(3) 



where it is implied that the lul and are the laser 
pulse carrier frequency and wave- vector, respectively and 
K' = dK/dhJ is related only to material dispersion and 
is evaluated at w^. In general we may thus indepen- 
dently choose the angles 7 and 5 so that the phase and 
peak velocities are effectively two independent and tun- 
able quantities. We focus our attention on the specific 
case 5 = —7, i.e. of the PBB, as this allows sufficient 
control over the phase and peak velocity difference whilst 
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FIG. 2: (in color online) XUV relative spectral energy (in log- 
arithmic scale) versus photon energy for three different PBB 
cone angles: dashed curve - 1.5 deg, dotted curve - 1 deg, solid 
curve - 0.5 deg. The thin black line indicated with deg is 
the reference spectrum obtained with a Gaussian pump pulse 
that has a FWHM of 40 nm. The spectral energy for 1 and 
1.5 deg has been normalized to that of 0.5 deg at the cutoff 
(116.2 eV). The vertical dashed lines indicate the positions of 
the 21*^* and 41"* harmonics, analyzed in more detail in Fig.|3] 



remaining experimentally simple to access. For the case 
of a PBB, Eq. ^ becomes Vp = cos-f/K'. The difference 
between the phase and peak velocities will induce a beat- 
ing in the pulse peak intensity with a periodicity that is 
twice that of the CEP and is a function of the cone angle 
7, 



A(7) = 



(4) 



So, for example, for a Bessel cone angle 7=1 deg, 
Eq. ([4]) predicts a periodicity in the pulse intensity 
modulation A = 1.28 mm. 

In Fig. [TJa) we show a cross-section of the geometry 
of the PBB. The tilted wavefronts propagate at an angle 
with respect to the propagation direction, z. The region 
in which they overlap is called the Bessel zone (darker 
shaded area in the figure) : the line corresponding to z=0 
indicates the position at which the PBB is considered in 
the present study as the input condition. Such a pulse 
and input condition may be obtained for example by 
illuminating a circular diffraction grating with a colli- 
mated Gaussian beam [Tsl . [l^ or by nonlinear frequency 
conversion processes in nonlinear optical crystals [20[ . 

In order to verify the validity of this novel QPM 
technique we performed numerical simulations based on 
a three dimensional propagation model in cylindrical 
coordinates using the nonadiabatic Strong-Ficld- 
Approxiniation l2lj to calculate the atomic response, 
as outlined in [6|. The computational code is fully 
parallelized with typical runs of 24 hours using 256 
CPU in a high performance interprocess communication 
system. The laser pulse central wavelength was chosen 
to be 800 nm, and the input field was written as E = 
EoJo{rojL sin7/c) exp(-(r/p)^) exp(-(t/r)2) cos(tJi,i) 
with a pulse duration at full-width-half-maximum 
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FIG. 3: (in color online) 21''' and harmonic relative spec- 
tral energy (in logarithmic scale) versus propagation distance 
for two different PBB cone angles, 7 = 1, 1.5 deg. 



(FWHM) of 5 fs (approximately 2 optical cycles). The 
Gaussian apodization in the transverse dimension was 
chosen to have a FWHM of 1.2 mm and the peak input 
intensity in all cases was 5 x 10^^ W/cm^. HHG was 
performed in Neon gas at a pressure of 20 mbar over a 
propagation distance of 5 mm. These values were chosen 
so as to have very weak reabsorption and thus observe 
the effect in the simplest case. In Fig. [TJb) we show 
the evolution along the propagation direction of the 
peak intensity of the PBB for two different cone angles. 
It is the CEP-induced oscillations visible in the figure 
that drive the QPM process with periodicities that 
correspond to Eq. (|4|). The slow decrease in the average 
intensity is due mainly to effective negative second order 
dispersion which arises due to the pulse front tilt in 
the PBB JJJ. This negative dispersion at the pump 
pulse frequency may be partially compensated for by the 
positive material dispersion or eliminated by resorting to 
the use of more complicated non-dispersive pulses known 
as X- waves [13, 112] ■ However, not withstanding this slow 
intensity decrease, clear enhancement of phase-matching 
is observed as shown in Fig. [2l 

Three different cases are considered, 7 = 0.5, 1 and 
1.5 deg. The radially integrated spectral energy of the 
XUV radiation (in Joules) per photon energy unit (in 
eV) is shown in Fig. [2]for 7 = 0.5, 1, 1.5 deg (the energies 
have been normalized to that of 7 = 0.5 deg at the 
cutoff (116.2 eV) for clarity). All three PBBs show a 
nearly 2 decade enhancement of the XUV intensity with 
respect to that of a Gaussian pulse. Most interestingly 
it is also clear how changing the PBB cone angle 
gradually shifts the phase-matched spectral region (with 
increasing photon energies for decreasing cone angles). 
Furthermore, the bandwidth of the quasi-phase-matched 
spectral region increases with decreasing cone angle. 

In Fig. [3] we show the XUV relative spectral 
energy evolution versus the propagation distance for 
two different harmonics, the 21*^' and the 41^' (whose 
spectral location is indicated by the dashed vertical 
lines in Fig. [2]) and for two different PBB cone angles, 
1 and 1.5 deg. These harmonics are chosen so that one 
is well within the phase-matched spectral region and 
the other is outside. As can be seen, the quasi-phase- 
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FIG. 4: 21^* (in color online) Black lines - harmonic energy 
for two different values of the initial CEP, (a) and 7r/2 (b). 
The red lines show the laser pulse intensity (right axis) for 
the two cases. The shaded areas indicate the periodic plateau 
regions in the XUV QPM process. 



matched 21''* harmonic intensity is clearly increasing 
over the whole propagation distance whilst the non 
quasi-phase-matched 41*** harmonic intensity continues 
to oscillate around a nearly constant value. This gives 
clear evidence of a phase-matching mechanism in the 
intensity-enhanced spectral regions. 

In Fig. m we examine the XUV signal propagation in 
more depth. In particular we study the dependence of 
the 21'^' harmonic on the CEP oscillations and on the 
input CEP value for a PBB with 7 = 1 deg. In Fig. H^a) 
the red line shows the laser pulse on-axis intensity for an 
initial CEP=0 and the black line shows the correspond- 
ing XUV energy evolution in linear scale. After an initial 
relatively monotonic growth, the XUV signal shows a 
periodic pattern in which the energy increases, reaches 
a plateau, increases once more, etc., hence clearly 
demonstrating the signature of QPM. The shaded areas 
highlight the step-like behavior in the XUV signal that 
occurs with the same periodicity, A 1.28 mm, of the 
CEP induced oscillations and correspond to well-defined 
recurring positions in the periodic intensity pattern. In 
order to further investigate the influence of the CEP on 
phase-matching we performed an identical simulation 
with a different input CEP at z=0, CEP=7r/2. This is 
shown in Fig. IDJb) . The resulting shift by half a period 
in the laser intensity modulation delays the onset of the 
XUV energy growth by a corresponding half period. The 
following XUV signal evolution is then nearly identical 
to that observed for CEP=0 albeit with the same 
half-period shift to longer z. The harmonic generation 
is therefore strongly determined by the CEP induced 
oscillations and is locked not only to its periodicity but 
also to the precise value at each z. Finally, we note that 
the intervals in which the XUV energy increases do not 
coincide with the laser-pulse intensity maxima. This 
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indicates that the periodic increase in the XUV signal 
cannot be simply ascribed to an increase in the laser 
intensity but is rather related to a periodic modulation 
in the phase contributions to the nonlinear process. It is 
precisely this periodic phase modulation that gives rise 
to the observed QPM in the harmonic energy. 

In conclusion we report a novel QPM mechanism by 
which HHG may be enhanced in a tunable XUV spectral 
region. This mechanism relies on the fact that few-cycle 
laser pulses with sufficiently different phase and group 
velocities will undergo periodic intensity oscillations due 
to CEP evolution during propagation. In the present 
work the phase and group velocities (and thus the CEP) 
arc controlled by shaping the laser pulse into a pulsed 
Bessel beam and by then varying the Bessel cone angle. 
The resulting intensity oscillations arc shown to give 
rise to QPM and a nearly two decade enhancement in 
the XUV intensity. This technique is very general and 



may be applied to other nonlinear frequency conversion 
regimes once a suitable method for controlling the CEP 
is adopted. 
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